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FOREWORD 

This  work  supports  an  ongoing  effort  at  the  Naval  Surface  Warfare  Center,  Dahlgren 
Division  (NSWCDD)  to  determine  effective  algorithms  for  correlating  tracks  from  radar  and 
Electronic  Support  Measures  (ESM)  systems  with  angular  measurement  accuracies  on  the  order 
of  0.1  deg.  This  research  investigated  two  tracking  filters  that  can  be  used  in  an  ESM  system 
to  generate  bearing  and  elevation  tracks  and  was  performed  in  support  of  the  Surface-Launched 
Weaponry  Technology  Block. 

This  report  was  reviewed  by  Larry  Sumner,  Electronic  Warfare  Systems  Integration 
Branch,  Tom  Kimbrell,  Head,  Electronic  Warfare  Systems  Integration  Branch;  and  Richard  Lee, 
Head,  Electronic  Warfare  Systems  Division. 
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ABSTRACT 


Both  an  a,Pand  an  a, P,y  filter  have  been  used  to  estimate  the  true  bearing  and  elevation 
of  high-speed  maneuvering  targets  from  measurements  made  by  a  shipboard  electronic  support 
measures  (ESM)  sensor.  The  motivation  for  this  research  is  the  need  for  correlating  radar  and 
ESM  tracks.  Four  Mach  3  target  trajectories  with  a  data  rate  of  10  Hz  and  a  measurement  noise 
standard  deviation  of  0.1  deg  were  used  to  test  the  filters.  Consistency  tests  were  used  to 
determine  the  best  values  for  the  steady-state  gains  of  the  filters.  Both  filters  are  computationally 
simple  and  provide  estimates  of  greater  accuracy  than  the  ESM  measurements.  It  was  determined 
that  an  a, p  filter  provides  better  estimates  than  an  a, p,y  filter  for  most  of  the  tested  target 
trajectories.  Two  exceptions  are  the  bearing  for  a  Mach  3  crossing  target  at  ranges  less  than  20 
km  and  the  elevation  of  a  Mach  3  target  at  ranges  less  than  5  km.  Both  filters  are  recommended 
for  radar-to-ESM  track  correlation  research.  The  a, p  filter  should  be  tested  first  since  it  is 
simpler  and  provides  better  estimates  for  most  of  the  cases  tested. 
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CHAPTER  1 
INTRODUCTION 

Electronic  support  measures  (ESM)  sensors  detect  and  make  measurements  on 
electromagnetic  energy  emitted  from  radar  systems  in  the  environment;  these  measurements 
include  frequency,  pulse  width,  pulse  amplitude,  pulse  repetition  interval,  and  intra-pulse 
modulations.  However,  the  only  measurements  that  can  be  used  to  track  the  kinematics  of  the 
target  are  bearing  and  elevation. 

The  primary  requirement  for  having  backing  filters  in  a  shipboard  ESM  system  comes 
from  the  need  to  correlate  ESM  tracks  with  the  tracks  of  radars  and  other  onboard  sensors.  It  is 
not  necessary  for  the  ESM  system  to  be  able  to  determine  the  target’s  position  in  three- 
dimensional  space  since  the  radar  tracks  provide  this  information.  However,  it  is  necessary  for 
the  ESM  system  to  provide  accurate  angle  rate  estimates  in  order  to  reliably  correlate  its  tracks 
with  those  of  the  radar,  so  that  these  systems  can  be  used  cooperatively  against  all  threats. 

Another  requirement  for  the  ESM  system  is  to  track  hundreds  of  platforms  at  once,  which 
will  force  the  tracking  filters  to  be  computationally  simple.  For  this  reason,  this  research 
investigated  two  of  the  more  computationally  simple  tracking  filters  known.  Further  research  will 
be  performed  to  determine  the  effectiveness  of  these  filters  for  correlating  radar  and  ESM  tracks. 
If  these  filters  prove  inadequate  for  this  function,  then  other  more  complex  tracking  filters  will 
be  investigated  for  this  purpose. 

A  single  type  of  tracking  filter  will  not  likely  provide  adequate  performance  for  all  targets 
since  an  ESM  sensor  must  track  many  targets  with  widely  varying  capabilities  for  maneuvering 
and  with  widely  varying  pulse  repetition  intervals.  However,  ESM  sensors  can  associate 
measurements  with  the  correct  target  with  a  high  probability  of  success  due  to  the  large  number 
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of  parameters  they  measure  and  their  target  identification  abilities.  Therefore,  the  target  tracking 
function  for  an  ESM  sensor  can  use  the  emitter  identification  information  to  assign  an  appropriate 
tracking  filter  for  each  target  based  on  the  platform  type;  e.g.,  the  ESM  tracking  algorithm  can 
have  a  table  of  filter  gains  and  data  rates  based  on  platform  type  (missile,  attack  aircraft,  ships, 
etc.).  Then,  the  appropriate  gains  and  data  rates  can  be  chosen  to  track  that  particular  platform. 

The  research  discussed  here  concentrates  on  evaluating  the  performance  of  two 
computationally  simple  tracking  filters  against  highly  maneuverable  Mach  3  targets  using  an 
extremeiy  slow  data  rate  (i.e.,  on  the  order  of  10  Hz).  For  slower  targets,  both  filters  can  be 
adjusted  to  provide  better  estimates  by  decreasing  the  process  noise  variance  (i.e.,  the  steady  state 
gains).  Likewise,  if  the  ESM  sensor  can  give  a  higher  data  rate,  then  the  tracking  filter  can  be 
adjusted  to  improve  the  quality  of  the  estimates.  Therefore,  further  research  must  be  performed 
using  slower  targets  and  faster  data  rates  in  order  to  obtain  optimal  performance  from  the 
tracking  filters  discussed  in  this  report. 

The  tracking  filters  for  radar  sensor  applications  have  generally  been  implemented  in 
Cartesian  coordinates.  This  is  possible  because  the  radar  measures  all  three  sphericai  coordinates 
(range,  bearing,  and  elevation),  so  a  transformation  to  Cartesian  coordinates  can  be  performed. 
Most  shipboard  ESM  sensors  do  not  make  range  measurements  because  to  do  so  requires  either 
a  maneuver  of  the  sensor  platform  or  triangulation  of  angle  measurements  from  two  or  more 
sensors.  This  research  assumed  that  no  range  measurement  is  available  since  the  targets  of 
interest  can  perform  maneuvers  of  greater  magnitude  than  a  ship.  As  discussed  by  Gorecki,  range 
is  unobservable  by  a  passive  sensor  unless  the  sensor  platform  performs  maneuvers  (in  directions 
orthogonal  to  the  line  of  sight)  of  at  least  one  order  of  magnitude  greater  than  the  maneuvers 
performed  by  the  target.1  Also,  triangulation  will  likely  require  a  data  link  between  ships  which 
may  not  be  available  at  all  times. 

A  tracking  algorithm  that  has  been  successfully  implemented  in  aircraft  passive  sensor 
applications  is  the  Modified  Spherical  Coordinates  (MSC)  angle-only  filter.2,3  It  uses  the  Extended 
Kalman  Filter  and  states  that  are  defined  in  the  MSC  coordinate  system  in  a  manner  that 
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decouples  the  range  estimate  from  the  angle  and  angle  rate  estimates  to  prevent  degradation  of 
the  quality  of  the  angle  and  angle  rate  estimates  when  good  range  estimates  are  not  available.  The 
aircraft  can  only  obtain  good  range  estimates  using  passive  sensors  by  maneuvering  when  the 
target  is  not  maneuvering  or  by  triangulation  using  measurements  from  the  passive  sensors  of 
other  platforms.  If  effective  and  reliable  techniques  for  obtaining  range  estimates  using  shipboard 
ESM  sensors  are  developed,  then  the  MSC  angle-only  filter  should  be  considered.  However,  the 
lack  of  passive  ranging  techniques  and  the  computational  complexity  of  the  MSC  angle-only  filter 
make  it  an  unlikely  candidate  to  accomplish  the  goals  of  this  research. 

An  a,p  and  an  a,P,y  filter  have  been  tested  as  ESM  tracking  filters  against  four  different 
target  trajectories.  Both  filters  track  bearing  and  elevation  measurements  independently.;  Both 
filter  designs  are  based  on  the  assumption  that  the  true  bearing  and  elevation  measurements  are 
uncorrelated.  The  aue  bearing  and  elevation  will  be  correlated  with  range;  but,  as  stated  above, 
it  is  assumed  that  no  range  measurements  are  available.  Consistency  tests  are  used  to  determine 
appropriate  values  for  the  process  noise  variances.  Results  are  presented  to  demonstrate  that  both 
filters  provide  adequate  bearing  and  elevation  estimates  for  all  of  the  tested  trajectories. 
Significance  tests  of  the  differences  in  mean  squared  errors  are  used  to  compare  the  performance 
of  the  a,P  and  the  a.P.y  filters.  It  is  shown  that  the  a,p  filter  provides  better  estimates  for  most 
of  the  trajectories  tested. 
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CHAPTER  2 
FILTER  DESCRIPTIONS 


The  two  tracking  filters  tested  are  both  standard  Kalman  filters  operating  on  the  bearing 
and  elevation  measurements  independently.  The  only  difference  in  the  two  tracking  filters  is  the 
number  of  states  estimated;  the  a, P  filter  uses  a  constant  first  derivative  prediction  matrix 
(referred  to  as  a  2  state  filter),  and  the  a,  p,y  filter  uses  l  jnstant  second  derivative  prediction 
matrix  (referred  to  as  a  3  state  filter).  The  Kalman  filter  equations  are  presented  first  followed 
by  the  prediction  matrices  for  the  2  and  3  state  filters. 


The  Kalman  filter  equations  are  presented  in  order  to  define  several  terms  that  will  be 
used  in  later  sections.  The  first  step  performed  by  a  Kalman  filter  is  to  predict  the  current 
estimate  and  its  covariance  to  the  time  of  the  next  measurement  using  the  equations 

(Z2> 

where  is  the  state  estimate  at  time  j  given  the  measurements  up  to  time  k-f*  is  its 

covariance,  Fk  is  the  state  transition  matrix,  and  Qk  is  the  covariance  of  the  process  noise.4  The 

next  step  is  to  calculate  the  covariance  of  the  innovations  and  then  the  filter  gain  using  the 
equations 


(2.3) 

(2.4) 


where  Hk  is  the  measurement  matrix,  Rk  is  the  covariance  of  the  measurement  noise,  Sk  is  the 
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covariance  of  the  innovation  (defined  below),  and  Wk  is  the  filter  gain.4  The  last  step  is  to 
calculate  the  innovation,  the  .iew  state  estimate,  and  its  covariance  using  the  equations 


V*+l  ~Zk*  1  ~^k*  A+l|fc 

(2.5) 

■4+l|*+l  ^*+lV*+l 

(2.6) 

*+lHr+l  =^k*  11*"  ^k*\ 

(2.7) 

where  zk  is  the  measurement  at  time  k  and  vk  is  the  innovation.4 


These  Kalman  filter  equations  are  used  for  both  the  2  and  3  state  tracking  filters.  The  only 
differences  between  these  filters  are  the  matrices  Fk  and  Qk  in  Equations  (2.1)  and  (2.2)  above. 

For  the  2  state  filter,  Fk  is  given  by5 


1  J\ 


P  iJ 


(2.8) 


where  T  is  the  time  step  between  measurements;  for  the  3  state  filter,  Fk  is  given  by5 


1  t  -T2 
2 

0  1  T 
0  0  1 


(2.9) 


The  process  noise  covariance  matrix,  Qk  ,  is  given  by 

Q^ohfik  OJ«) 

where  qk  is  the  process  noise  variance  that  must  be  set  by  the  filter  designer,  and  the  vector, 
Gk,  corresponds  to  the  Fk  for  the  filter.5  For  the  2  state  filter  with  Fk  in  Equation  (2.8),  Gk 
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is  given  by 

G.Jlr2  jj 

and  by5 

Gk=-T 2  T  1 

[2 

for  the  3  state  filter  with  Fk  in  Equation  (2.9).5 


(2.11) 


(2.12) 
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CHAPTER  3 
CONSISTENCY  TESTS 


Several  consistency  tests  can  be  used  to  match  the  process  noise  variance,  qk,  in  Equation 

(2,10)  to  a  particular  target  trajectory.  The  consistency  tests  that  are  performed  for  Monte  Carlo 
simulations  are  described  first,  followed  by  descriptions  of  the  tests  that  can  be  performed  on  a 
single  run. 


The  consistency  tests  for  Monte  Carlo  simulations  of  tracking  filters  are  performed  on  the 
sample  average  of  the  normalized  state  estimation  error  squared,  the  sample  average  of  the 
normalized  innovation  squared,  and  the  sample  autocorrelation  statistic  of  the  innovations.  The 
normalized  state  estimation  error  squared  is  given  by 

(3,1) 

where 


(3.2) 


is  the  difference  between  the  true  state  and  the  state  estimate  at  time  k  4  The  sample  average 
of  is  computed  over  all  of  the  Monte  Carlo  runs  by  the  equation4 


i  N 
iV  /.  i 


(3.3) 


The  sample  average,  et,  is  compared  to  a  confidence  region  defined  by  thresholds  obtained  from 
a  Chi-squared  table.  If  ek  is  within  the  confidence  region,  the  filter’s  errors  are  consistent  with 
the  filter’s  covariances.  A  similar  test  is  performed  on  the  sample  average  of  the  normalized 
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innovation  squared.  The  normalized  innovation  squared  is  computed  by  the  equation4 


4=v*Vv* 


(3.4) 


The  sample  average  of  evk  is  denoted  by  e£.  The  last  consistency  test  for  Monte  Carlo  simulations 
uses  a  sample  autocorrelation  statistic  of  the  innovations  that  is  computed  by  the  equation4 


N 


ivE&W 

i»  1 


N 


N 


E 

i«X  i-1 


(3.5) 


The  filter’s  innovations  are  considered  uncorrelated  (or  white)  if  Pjy  is  within  a  confidence 
region  defined  by  thresholds  taken  from  a  table  of  the  Normal  density  function. 


The  consistency  tests  of  the  innovations  can  also  be  performed  on  a  single  run  of  the  filter 
by  using  time  averages  instead  of  sample  averages.  The  test  statistic  for  consistent  innovations 
in  a  single  run  is 


^=iEv*Vv* 


(3.6) 


and  the  test  statistic  for  the  autocorrelation  of  the  innovations  in  a  single  run  is4 


k-l 


p,=E  (v/vJE  (v/v*E 


k- 1 


*-l 


(3.7) 
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CHAPTER  4 
FILTER  COMPARISONS 


A  comparison  between  two  tracking  Filters  can  be  made  using  a  significance  test  of  the 
sample  mean  of  the  differences  in  squared  error.  Comparing  the  root-mean-squared  errors 
(RMSE)  alone  is  not  effective  because  a  confidence  region  is  needed  to  determine  if  the 
differences  in  RMSE  are  significant.  The  hypothesis  that  algorithm  1  is  better  than  algorithm  2 
is  accepted  if  the  sample  mean  of  the  differences  in  squared  error  given  by 


(4.1) 


is  greater  than  a  threshold  taken  from  a  table  of  the  Normal  density.4  In  Equation  (4.1), 


(4.2) 


is  the  sample  mean  of  the  differences  in  squared  error  where  C"  is  the  squared  error  for 
algorithm  m  at  a  particular  time  and 


is  the  variance  of  the  differences  in  squared  error.4 


(4.3) 
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CHAPTER  5 

FILTER  INITIALIZATION 


The  procedure  for  choosing  the  initial  estimates  and  their  covariances  is  described  in  this 
section.  For  the  2  state  filter,  the  initial  estimate  is  computed  from4 

V[Zo  fco -*-,V7]r  <5-'> 

Similarly,  the  initial  estimate  for  the  3  state  filter  is  computed  from 


■*0(0  I 


*0 

(Zq-z.JIT 

[(zo-Zz-.n.j/r1 


The  covariances  are  given  by4 


P  = 
^ojo 


*0  V 

RJT  2 RJT2 


for  the  2  state  filter  and  by 


0K> 


*o  RJT  RJT 2 
RJT  2 RJT2  3RJT* 
RJT2  3RJT3  6RJT 4 


for  the  3  state  filter. 


(5.2) 


(5.3) 


(5.4) 
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CHAPTER  6 

TARGET  TRAJECTORIES 

ESM  measurements  can  be  generated  by  converting  trajectories  formed  in  Cartesian 
coordinates  into  spherical  coordinates  and  adding  pseudorandom  noise  sequences  to  the  bearing 
and  elevation  data.  Two  unrealistic  but  reasonable  assumptions  make  the  generation  of  the  target 
trajectories  easier.  One  is  to  assume  the  earth  is  flat,  which  is  reasonable  for  the  short  ranges 
used  in  the  trajectories  discussed  in  this  report;  the  other  is  to  assume  the  ESM  sensor  is  at  zero 
altitude.  This  assumption  causes  the  true  elevation  measurements  to  always  be  positive.  For  an 
actual  shipboard  system  tracking  a  sea  skimming  missile,  the  true  elevation  measurements  will 
be  negative.  However,  the  magnitude  of  the  true  elevation  rate  is  the  same  in  either  case,  so  the 
positive  elevation  measurements  are  considered  to  be  reasonable.  Zero  mean  pseudorandom 
number  sequences  with  a  white  Gaussian  distribution  are  used  to  add  errors  to  the  true  bearing 
and  elevation.  A  standard  deviation  of  0.1  deg  (1.73  mrad)  and  a  data  rate  of  10  Hz  is  used  for 
all  the  simulations  discussed  here.  For  Monte  Carlo  simulations,  a  new  pseudorandom  number 
sequence  is  generated  for  each  run,  so  the  measurement  errors  are  independent  from  run  to  run. 

Four  different  target  trajectories  were  generated  to  test  the  ESM  tracking  filters;  all  are 
Mach  3  targets  that  remain  at  constant  altitude.  The  ESM  sensor  is  located  at  the  origin  for  all 
four  trajectories.  The  first  target  moves  in  a  straight  line  directly  towards  the  sensor;  it  is  referrea 
to  as  the  radial  target.  The  second  target  moves  in  a  straight  line  on  a  crossing  trajectory;  it  is 
referred  to  as  the  crossing  target.  The  third  target  moves  in  a  straight  line  towards  the  sensor  and 
performs  a  single  S-tum  in  the  middle  of  the  trajectory;  it  is  referred  to  as  the  S-turn  target.  The 
last  target  is  a  radial  target  that  performs  jinking  type  maneuvers;  it  is  referred  to  as  the  jinking 
target.  Each  of  these  trajectories  is  presented  (all  of  the  figures  referred  to  in  this  section  are 
located  in  Appendix  A). 
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6.1  RADIAL  TARGET 

The  radial  target  starts  at  XY  coordinates  (25  km,  0  m)  and  ends  at  (510  m,  0  m).  The 
true  bearing,  elevation,  and  elevation  rate  for  this  target  are  shown  in  Figures  A-l,  A-2,  and  A-3, 
respectively.  The  target  maintained  a  constant  altitude  of  20  m  and  a  constant  speed  of  Mach  3 
(1029  m/s)  throughout  the  trajectory.  The  elevation  increases  almost  linearly  until  about  20  s  into 
the  trajectory,  where  it  increases  rapidly  because  the  target  is  at  close  range. 

6  2  CROSSING  TARGET 

A  plot  of  the  XY  plane  view  of  the  crossing  target  is  given  in  Figure  A-4.  The  crossing 
target  starts  at  XY  coordinates  (10  km,  -23  km)  and  ends  at  (10  km,  23  km).  These  points  are 
labeled  *t=0’  and  ‘1=41.0’  in  Figure  A-4,  The  closest  point  of  approach  occurs  midway  through 
the  trajectory  at  a  range  of  10  km.  The  target  maintained  a  constant  altitude  of  20  m  and  a 
constant  speed  of  Mach  3  (1029  m/s)  throughout  the  trajectory.  The  true  bearing  and  bearing  rate 
for  this  target  are  shown  in  Figures  A-5  and  A-6,  respectively.  As  shown  is  Figure  A-6,  the 
largest  absolute  value  of  the  bearing  rate  occurs  at  20  s  into  the  trajectory,  and  the  maximum 
change  in  the  bearing  rate  occurs  at  10  and  30  s  into  the  trajectory.  The  true  elevation  and 
elevation  rate  for  this  target  are  shown  in  Figures  A-7  and  A-8,  respectively. 

6.3  S-TURN  TARGET 

A  plot  of  the  XY  plane  view  of  this  trajectory  is  given  in  Figure  A-9,  The  target  starts 
at  XY  coordinates  (25  km,  -925  m),  and  ends  at  (1.95  km,  -58  m).  The  target  maintained  a 
constant  altitude  of  30.5  m  and  a  constant  speed  of  Mach  3  (1006  m/s)  throughout  the  trajectory. 
The  target  moves  in  a  straight  line  for  10  s.  It  then  performs  a  15  g,  constant  speed  turn  towards 
the  positive  Y  direction  for  2.5  s  and  then  a  15  g,  constant  speed  turn  towards  the  negative  Y 
direction  for  2.5  s.  The  target  moves  in  a  straight  line  towards  the  ESM  sensor  for  the  last  8  s 
of  the  trajectory.  The  velocities  and  accelerations  for  this  trajectory  are  shown  in  Figures  A- 10 
and  A-ll.  Notice  in  Figure  A- 11  that  there  are  no  instantaneous  changes  in  the  acceleration 
because  they  were  filtered  through  a  second-order  system.  The  true  bearing,  bearing  rate, 
elevation,  and  elevation  rate  are  shown  in  Figures  A-12,  A-13,  A-14,  and  A-15,  respectively. 
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6.4  JINKING  TARGET 

A  plot  of  the  XY  plane  view  of  this  trajectory  is  given  in  Figure  A- 16.  The  target  starts 
at  XY  coordinates  (25  km,  0  m),  and  ends  at  (1.06  k  n,  0  m).  These  points  are  labeled  *t=0.0’ 
and  ‘t=23.8’  in  Figure  A-16.  Also,  every  10th  data  point  is  labeled  with  an  ‘O’  in  the  figure  for 
clarity.  The  target  maintained  a  constant  altitude  of  50  m  and  a  constant  speed  of  Mach  3  (1006 
m/s)  throughout  the  trajectory.  The  target  moves  in  a  straight  line  for  5  s.  It  then  begins  jinking 
maneuvers  that  last  until  21.8  s  through  the  trajectory.  The  term  jinking  maneuver  refers  to 
alternating  accelerations  orthogonal  to  the  velocity  vector.  The  velocities  and  accelerations  for 
this  target  are  shown  in  Figures  A- 17  and  A- 18,  respectively.  The  Z  velocity  for  this  target  is 
zero  and  is  not  shown  in  Figure  A- 17.  Likewise,  the  X  and  Z  accelerations  are  zero  and  are  not 
shown  in  Figure  A- 18.  The  Y  acceleration  changes  instantaneously  from  +15  g  to  -15  g.  The 
accelerations  were  not  filtered  through  the  second-order  system  to  make  it  easier  to  direct  the 
target  toward  the  ESM  sensor  located  at  the  origin.  The  target  moves  in  a  straight  line  towards 
the  ESM  sensor  for  the  last  2  s  of  the  trajectory.  The  true  bearing,  bearing  rate,  elevation,  and 
elevation  rate  are  shown  in  Figures  A- 19,  A- 20,  A-21,  and  A-22,  respectively.  Notice  in  Figures 
A- 19  and  A-20  that  the  jinking  maneuvers  cause  faster  changes  in  the  bearing  as  the  target  closes 
in  range. 
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CHAPTER  7 

PROCESS  NOISE  VARIANCE  SELECTION 

The  filter  consistency  tests  are  used  to  select  the  best  value  of  the  process  noise  variance 
for  each  filter  for  each  trajectory.  The  single  run  tests  are  used  to  determine  a  reasonable  range 
of  values  for  the  process  noise  variance,  and  the  Monte  Carlo  simulations  are  performed  to  find 
the  best  value  for  the  variance.  An  explanation  of  the  procedures  used  to  select  the  thresholds 
for  the  consistency  tests  are  given.  Then,  the  variances  chosen  for  each  filter  for  each  trajectory 
are  presented.  This  chapter  concludes  with  an  example  of  the  procedure  used  to  select  the  process 
noise  variance  for  a  filter. 

7.1  CONSISTENCY  TEST  THRESHOLDS 

The  sample  average  of  the  normalized  state  estimation  error  squared,  ek,  given  by 
Equation  (3.3),  is  tested  during  Monte  Carlo  simulations.  This  statistic  is  compared  to  a  threshold 
from  a  Chi-square  table  and  has  Nnx  degrees-of-freedom  (DOF)  where  N  is  the  number  of 

Monte  Carlo  runs  and  nx  is  the  number  of  states  estimated  by  the  filter.4  The  95-percent 
confidence  region  for  the  2  state  filter  (a  95-percent  confidence  region  was  used  in  the  tests  for 
all  statistics)  is  1.5  to  2.6  since  there  are  100  DOF  (N=5G  and  nx= 2)  and  the  limits  from  the 

table  (74.2  and  129.6)  must  be  divided  by  N.  Similarly,  the  confidence  region  for  the  3  slate 
filter  is  2.2  to  3.9  since  there  are  90  DOF  (N= 30  and  ^=3)  and  the  limits  from  the  table  (65.6 
and  118.1)  must  be  divided  by  N. 
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The  normalized  innovations,  e*  given  by  Equation  (3.4),  is  tested  during  Monte  Carlo 
simulations  and  single  runs  of  the  filters.  Both  test  statistics  are  Chi-squared  distributed,  but  with 
a  different  number  of  DOF.  The  sample  average  statistic,  Tk,  has  Nnz  DOF  where  nz  is  the 

number  of  measurements  for  each  estimate  (1  for  all  cases  considered  here).  The  confidence 
region  for  the  2  state  filter  is  0.65  to  1.43  since  there  are  50  DOF  (N=5 0  and  nz= 1 )  and  the 
iimits  from  the  table  are  32.3  and  71.4.  Similarly,  the  confidence  region  for  the  3  state  filter  is 
0.56  to  1.57  since  there  are  30  DOF  (N= 30  and  «z=l)  and  the  limits  from  the  table  are  16.8 

and  47.0.  The  time  average  statistic,  ev  given  by  Equation  (3.6),  has  Knz  DOF,  where  K  is  the 

number  of  estimates  used  to  compute  the  average.  The  confidence  region  for  both  the  2  and  3 
state  filters  is  0.74  to  1.3  since  the  time  averages  were  computed  using  blocks  of  100  estimates. 
Only  100  estimates  were  used  since  the  largest  degree  of  freedom  listed  in  the  Chi-squared  tables 
is  100.  For  DOF  greater  than  100,  the  statistic  is  approximately  Gaussian.4  However,  the  author 
chose  to  limit  the  number  of  estimates  used  to  compute  the  average  to  100  instead  of 
approximating  the  distribution  by  a  Gaussian. 

The  sample  autocorrelation  statistic,  given  by  Equation  (3.5),  is  tested  during  Monte 
Carlo  simulations,  and  the  time  average  autocorrelation  statistic,  p,  given  by  Equation  (3.7),  is 
tested  during  single  run  tests.  Onl}  the  single  step  difference  autocorrelation  statistic  is  computed 
in  all  cases  (i.e.,  k-j- 1  for  p^  ).  For  large  N,  is  approximately  normal  with  variance  1  IN, 
so  the  confidence  region  is  taken  from  a  table  of  the  Normal  disMbution.4  The  confidence  region 
for  pjy  is  ±  0.28  for  the  2  state  filter  and  ±  0.36  for  the  3  state  filter.  Likewise,  for  large  K,  pl 

is  approximately  normal  with  variance  1  /K.A  The  confidence  region  for  p;  is  ± 0.2  for  both  the 
2  state  filter  and  the  3  state  filter  since  both  use  averages  of  100  estimates. 


! 
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7.2  CHOSEN  VALUES  FOR  PROCESS  NOISE  VARIANCES 

The  value  of  the  process  noise  variance,  qk,  was  chosen  as  the  maximum  value  for  which 

the  Filter  did  not  violate  any  of  the  Monte  Carlo  simulation  consistency  tests.  This  is  the  most 
effective  method  for  selecting  the  process  noise  variance  since  the  process  noise  accounts  for  all 
of  the  modeling  inaccuracies  of  the  Kalman  filter.  The  Kalman  filter’s  constant  state  transition 
matrix  is  the  dominate  cause  of  these  modeling  inaccuracies.  The  test  of  ek  was  relaxed  to 

include  values  smaller  than  the  lower  threshold  defining  the  confidence  region  since  this  was  the 
case  f or  both  filters  on  all  of  the  trajectories.  This  condition  indicates  that  the  filter’s  covariances 
predict  larger  errors  than  reality  (i.e.,  the  filter  is  pessimistic).  The  qk  chosen  for  each  filter  for 

each  trajectory  are  presented  in  Table  1.  The  units  for  all  the  values  shown  in  the  table  are 
rad2/s4. 


TABLE  1.  CHOSEN  VALUES  OF  PROCESS  NOISE  VARIANCE 


Trajectory 

2  State  Filte. 

3  State  Filter 

Bearing 

Elevation 

Bearing 

Elevation 

Radial 

103 

IO'3 

io-3 

10'3 

Crossing 

104 

104 

10'6 

IO'6 

S-Tum 

10‘3 

IO'7 

10'5 

IO'7 

Jinking 

5xl0'3 

5xl04 

T 

O 

X 

5x1  O'6 

7.3  EXAMPLE  OF  PROCESS  NOISE  VARIANCE  SELECTION 

The  and  p^-  of  the  bearing  for  30  Monte  Carlo  runs  of  the  3  state  filter  on  the 

crossing  trajectory  are  shown  in  Figures  B-l  and  B-2  which  are  located  in  Appendix  B.  In  these 
figures,  the  area  between  the  horizontal  lines  is  the  95-percent  confidence  region.  The  Y  axis 

label  in  Figure  B-l,  ‘SAMPLE  AVG  NORM  INNOV  -  BEAR’  denotes  the  e^for  the  bearing 
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innovations.  Likewise,  the  Y  axis  label  'SAMPLE  AUTOCORR INNOV  -  BEAR’  in  Figure  B-2 
denotes  the  p^for  the  bearing  innovations.  Notice  in  these  Figures  that  most  of  the  i*  and 

are  within  the  confidence  region  and  the  mean  of  Tk  over  the  trajectory  is  close  to  one  and  the 
mean  of  is  close  to  zero.  The  qk  for  this  filter  was  10'6.  Other  Monte  Carlo  simulations  for 

this  filter  were  performed  with  qk  of  10‘5  and  10'7.  For  <7*= IQ'5,  more  of  the  Tk  and  p^  were 
below  the  lower  threshold  of  the  confidence  region,  and  their  means  were  not  as  close  to  their 
optimum  values  of  1  and  0,  respectively.  For  ^=10"7,  a  significant  number  of  the  Tk  and  p^ 
exceeded  the  upper  threshold  of  the  confidence  region  at  the  middle  of  the  trajectory  where  the 
change  in  bearing  approaches  zero  as  shown  in  Figures  B-3  and  B-4.  Thus,  the  designer  selected  ^=10^ 
for  this  filter  and  trajectory.  This  procedure  was  repeated  for  all  four  trajectories  for  both  filters 
to  select  the  qk  shown  in  Table  1. 
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CHAPTER  8 
RESULTS 

The  techniques  described  were  utilized  to  evaluate  the  performance  of  2  and  3  state  filters 
for  four  target  trajectories.  The  results  of  the  2  state  filter  are  given  first  followed  by  the  results 
for  the  3  state  filter.  Finally,  comparisons  of  the  performance  of  the  2  and  3  state  filters  are 
presented. 

8.1  RESULTS  FOR  2  STATE  FILTER 

The  following  subsections  are  organized  in  the  same  manner.  Sample  bearing,  bearing 
rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the  2  state  filters  are  presented. 
Then,  the  RMSE  for  these  estimates  from  50  Monte  Carlo  runs  of  the  2  state  filters  are  presented. 
Each  subsection  concludes  with  plots  of  the  consistency  test  statistics  from  50  Monte  Carlo  runs 
of  the  2  state  filters.  All  of  the  figures  referred  to  in  this  section  are  located  in  Appendix  C. 

8.1.1  Radial  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
2  state  filters  on  the  radial  target  are  shown  in  Figures  C-l,  C-2,  C-3,  and  C-4,  respectively.  The 
true  states  are  plotted  in  each  figure  for  reference.  These  figures  are  intended  to  give  the  reader 
an  intuitive  understanding  of  the  quality  of  the  estimates  generated  by  the  filters.  Since  these 
estimates  are  from  a  single  run  of  the  filters,  no  conclusions  on  the  average  performance  of  the 
filters  should  be  drawn  from  them. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  50  Monte  Carlo 
runs  of  the  2  state  filters  on  the  radial  target  are  shown  in  Figures  C-5,  C-6,  C-l,  and  C-8, 
respectively.  The  horizontal  lines  in  Figures  C-5  and  C-l  correspond  to  the  standard  deviation 
of  the  measurement  noise  (0.1  deg  or  1.74  mrad)  and  is  provided  as  a  reference  to  determine  the 
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amount  of  noise  removed  by  the  filter.  In  Figure  C-5,  the  RMS  bearing  error  has  been  reduced 
from  1.74  mrad  to  1  mrad.  In  Figure  C-6,  the  RMS  bearing  rate  error  is  2.8  mrad/s.  In  Figure 
C-7,  the  RMS  elevation  error  has  been  reduced  from  1.74  mrad  to  1  mrad  through  most  of  the 
trajectory,  but  the  error  increases  rapidly  in  the  last  second  of  the  trajectory.  This  occurs  because 
the  filter  assumes  the  elevation  increases  linearly,  but  the  true  elevation  is  increasing  much  faster 
than  a  linear  slope.  The  same  effect  can  be  seen  in  the  elevation  rate  estimates.  The  author  is 
not  concerned  about  the  large  errors  in  the  elevation  and  elevation  rate  estimates  during  the  last 
second  of  the  trajectory  since  the  target  is  within  2  km  of  the  ship.  When  the  missile  reaches  a 
range  of  2  km  from  the  ship,  the  radar-to-ESM  correlation  should  have  already  been  completed. 
Also,  the  accuracy  of  the  elevation  estimates  will  still  be  sufficient  to  support  the  electronic 
countermeasures  (ECM)  system  since  it  will  use  beamwidths  that  are  much  larger  than  the 
elevation  errors. 

Plots  of  the  consistency  test  statistics  from  50  Monte  Carlo  runs  of  the  2  state  filters  on 
the  radial  target  are  given  in  Figures  C-9  through  C-14.  The  area  between  the  horizontal  lines  is 
the  95-percent  confidence  region.  The  Y  axis  labels  used  in  these  figures  are  defined  as  follows: 

•SAMPLE  AVG  NORM  INNOV’  denotes  i*  given  by  Equation  (3.4),  'SAMPLE  AUTOCORR 

INNOV’  denotes  given  by  Equation  (3.5),  and  'NORM  STATE  EST  ERROR  SQUARED’ 

denotes  given  by  Equation  (3.3).  The  Chi-square  test  of  the  bearing  innovations  is  given  in 

Figure  C-9.  Most  of  the  points  are  within  the  95-percent  confidence  region,  so  the  bearing 
innovations  are  consistent  with  their  covariances.  The  Chi-square  test  of  the  elevation  innovations 
is  given  in  Figure  C-10.  This  figure  shows  that  the  elevation  innovations  are  consistent  until  the 
last  seconds  of  the  trajectory.  For  reasons  mentioned  previously,  this  is  not  a  reason  for  concern. 
The  whiteness  test  of  the  bearing  innovations  is  shown  in  Figure  C-l  1.  <’or  this  figure  and  all 
subsequent  whiteness  tests,  only  the  for  k-j= 1  is  shown.  As  shown  in  this  figure,  the 

innovations  are  roughly  uncorrelated.  The  whiteness  test  for  the  elevation  innovations  is  shown 
in  Figure  C-12.  The  innovations  become  significantly  correlated  at  the  end  of  the  trajectory  due 
to  the  large  errors  in  the  elevation  estimates.  The  Chi-square  test  of  the  bearing  estimation  error 
is  shown  in  Figure  C-l 3.  The  confidence  region  is  not  shown  since  the  thresholds  are  on  the 
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order  of  one.  Since  the  test  data  are  much  smaller  than  the  lower  threshold,  the  author  concludes 
that  the  filter’s  error  covariances,  Pk  given  by  Equation  (2.2),  predict  larger  errors  than  actually 
occur  (i.e.,  the  filter  is  pessimistic).  This  condition  is  considered  acceptable  since  it  is  not  critical 
that  the  Pk  predict  the  exact  amount  of  error  in  the  estimates;  it  is  only  important  that  the  Pk 

do  not  underestimate  the  amount  of  error  in  the  estimates.  This  is  equivalent  to  using  a  one-sided 
confidence  region  instead  of  a  two-sided  confidence  region.  A  one-sided  confidence  region  will 
used  on  all  subsequent  Chi-square  tests  of  the  state  estimation  error.  The  Chi-square  test  of  the 
elevation  estimation  error  is  shown  in  Figure  C-14.  The  data  up  to  22  s  through  the  trajectory 
are  on  the  order  of  10'5  and  cannot  be  seen  in  the  figure.  The  large  errors  in  the  elevation  and 
elevation  rate  estimates  cause  the  sharp  increase  in  Hk  in  the  last  seconds  of  Figure  C-14.  Again, 

the  missile  is  at  a  very  short  range  at  this  point,  so  large  elevation  estimate  errors  are  not 
significant.  Thus,  all  of  the  consistency  tests  passed,  so  the  process  noise  variances  used  for  these 
simulations  are  good  values. 

8.1.2  Crossing  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
2  state  filters  on  the  crossing  target  are  shown  in  Figures  C-15  through  C-18.  In  Figure  C-15,  the 
difference  between  the  sample  bearing  estimate  and  the  true  bearing  is  too  small  to  be  seen.  The 
errors  in  the  sample  bearing  rate  estimate  in  Figure  C-16  are  smallest  in  the  middle  of  the 
trajectory  where  the  bearing  rate  is  nearly  constant  for  a  short  time.  This  makes  sense  because 
the  state  transition  matrix,  Fk  given  by  Equation  (2.8),  assumes  that  the  true  bearing  rate  is 

constant.  In  Figure  C-17,  the  sample  elevation  estimate  looks  similar  to  the  sample  bearing 
estimate  for  the  radial  target  shown  in  Figure  C-l  since  the  true  elevation  remains  almost 
constant  over  the  trajectory.  Similarly,  the  sample  elevation  rate  estimate  for  the  crossing  target 
shown  in  Figure  C-18  looks  almost  the  same  as  the  sample  bearing  rate  estimate  for  the  radial 
target  shown  in  Figure  C-2. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  50  Monte  Carlo 
runs  of  the  2  state  filters  on  the  crossing  target  are  shown  in  Figures  C-19,  C-20,  C-21,  and  C-22, 
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respectively.  In  Figure  C-19,  the  correlation  between  the  accuracy  of  the  state  transition  matrix, 
Fk ,  and  the  accuracy  of  the  estimates  is  demonstrated  clearly.  At  the  beginning,  middle,  and  end 

of  the  trajectory  where  the  true  bearing  rate  is  nearly  constant,  the  RMS  bearing  error  is  at  its 
minimum  of  0.8  mrad  (a  54-percent  reduction  from  the  1.74  mrad  measurement  noise  standard 
deviation).  At  the  15-s  and  27-s  points  of  the  trajectory,  the  bearing  rate  is  changing  at  its 
maximum  rate,  so  the  RMS  bearing  error  is  at  its  maximum  of  1.4  mrad  (only  a  19.5-percent 
reduction  from  1.74  mrad).  The  same  correlation  can  be  easily  seen  in  the  RMS  bearing  rate 
errors  shown  in  Figure  C-20.  The  minimum  RMS  bearing  rate  error  occurs  at  the  midpoint  of  the 
trajectory  at  1.25  mrad/s,  and  the  maximum  error  occurs  at  the  15-s  and  27-s  points  at  4.25 
mrad/s.  In  Figure  C-21,  the  RMS  elevation  error  has  been  reduced  from  1.74  mrad  to  0.8  mrad. 
The  RMS  elevation  rate  error  remains  at  1.25  mrad/s  throughout  the  trajectory  as  shown  in  Figure 
C-22. 


Plots  of  the  consistency  test  statistics  from  50  Monte  Carlo  runs  of  the  2  state  filters  on 
the  crossing  target  are  given  in  Figures  C-23  through  C-28.  The  Chi-square  test  of  the  bearing 
innovations  is  given  in  Figure  C-23.  Notice  that  a  large  number  of  the  points  are  outside  the  95- 
percent  confidence  region  at  the  15-s  and  27-s  points  of  the  trajectory  which  indicates  that  the 
bearing  innovations  are  not  consistent  with  their  covariances.  However,  if  the  process  noise 
variance,  qk,  is  increased  to  improve  the  consistency  of  the  innovations  at  these  two  points  of 

the  trajectory,  then  the  innovations  near  the  beginning,  middle,  and  end  of  the  trajectory  will 
become  less  consistent.  It  was  determined  through  many  Monte  Carlo  simulations  that  the  bearing 
innovations  for  this  trajectory  cannot  be  made  more  consistent  than  that  shown  in  Figure  C-23. 
In  Figure  C-24,  it  is  shown  that  the  elevation  innovations  are  very  consistent  through  the  entire 
trajectory.  This  is  not  surprising  since  the  true  elevation  rate  is  virtually  constant,  and  theFt 

predicts  a  constant  elevation  rate.  The  whiteness  test  of  the  bearing  innovations  is  shown  in 
Figure  C-25.  This  test  has  the  same  problems  that  were  encountered  for  the  Chi-square  test  of 
the  bearing  innovations  shown  in  Figure  C-23.  Again,  it  was  found  through  Monte  Carlo 
simulations  that  any  attempt  to  reduce  the  positive  correlation  of  the  bearing  innovations  by 
increasing  qk  will  cause  an  increase  in  the  negative  correlation  in  the  other  parts  of  the 
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trajectory.  The  whiteness  test  for  the  elevation  innovations  in  Figure  C-26  shows  that  they  are 
indeed  uncorrelated  in  time.  The  Chi-square  test  of  the  bearing  estimation  error  is  shown  in 
Figure  C-27.  The  bearing  estimation  error  is  consistent  using  a  one-sided  confidence  region. 
Using  a  two-sided  confidence  region,  this  filter  is  clearly  pessimistic.  The  Chi-square  test  of  the 
elevation  estimation  error  in  Figure  C-28  demonstrates  that  they  are  consistent.  As  in  the  radial 
target  case,  all  of  the  consistency  tests  passed,  so  the  process  noise  variances  used  for  these 
simulations  are  good  values. 

8.1.3  S-tum  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
2  state  filters  on  the  S-tum  target  are  shown  in  Figures  C-29  through  C-32.  In  Figure  C-29,  the 
reader  can  see  that  the  bearing  estimate  follows  the  true  bearing  maneuver  fairly  well,  but  the 
estimates  are  somewhat  noisy  before  and  after  the  S-tum  maneuver.  This  is  a  result  of  choosing 
the  qk  large  to  allow  significant  deviations  in  the  true  bearing.  The  same  is  true  of  the  bearing 

rate  estimate  in  Figure  C-30.  In  Figure  C-31,  the  sample  elevation  estimate  looks  similar  to  the 
sample  elevation  estimate  for  the  radial  target  shown  in  Figure  C-3.  However,  the  elevation 
estimate  in  Figure  C-31  lags  the  true  elevation  to  a  larger  degree  than  does  the  one  in  Figure  C-3. 
This  is  due  to  the  large  differences  in  the  qk  for  the  two  filters.  Referring  to  Table  1,  the  qk  for 

the  elevation  of  the  radial  target  was  chosen  as  10'3,  which  allows  the  filter  to  follow  the  changes 
in  the  true  elevation  more  closely  at  the  expense  of  degraded  estimates  when  the  true  elevation 
does  not  change.  The  qk  for  the  S-tum  target  was  chosen  as  10'7,  which  tightens  the  filter’s 

estimates  causing  it  to  lag  behind  the  changes  in  the  true  elevation.  Since  the  true  elevation  for 
these  two  targets  is  similar,  comparing  the  sample  elevation  estimates  from  these  two  filters 
provides  a  good  intuitive  understanding  of  the  relationship  between  qk  and  the  performance  of 
the  filter.  As  shown  in  Figure  C-32,  the  lag  in  the  elevation  rate  estimates  is  even  more  dramatic. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  50  Monte  Carlo 
runs  of  the  2  state  filters  on  the  S-tum  target  are  shown  in  Figures  C-33  through  C-36.  In  Figure 
C-33,  the  RMS  bearing  error  remains  nearly  constant  throughout  the  trajectory  at  about  1.1  mrad 
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(a  37-percent  reduction  from  the  1.74  mrad  measurement  noise  standard  deviation).  Notice  that 
one  cannot  easily  see  where  the  maneuver  occurs  near  the  middle  of  the  trajectory.  As  discussed 
in  the  previous  paragraph,  this  is  caused  by  the  rather  large  qk  chosen  for  this  filter  and 

trajectory.  In  Figure  C-34,  the  maneuver  is  easier  to  see  in  the  RMS  bearing  rate  errors.  If  th eqk 

where  smaller,  then  the  RMS  bearing  rate  error  before  and  after  the  maneuver  would  be  smaller, 
but  it  would  be  larger  during  the  maneuver.  In  Figure  C-35,  the  RMS  elevation  error  is  about  0.5 
mrad  in  the  middle  of  the  trajectory  which  is  a  reduction  of  71 -percent  from  1.74  mrad.  However, 
at  the  end  of  the  trajectory  the  elevation  errors  are  growing  quickly  which  is  the  consequence  of 
having  a  small  qk.  The  RMS  elevation  rate  error  (Figure  C-36)  is  very  small  in  the  middle  of  the 
trajectory  but  grows  quickly  at  the  end  as  did  the  RMS  elevation  error. 

Plots  of  the  consistency  test  statistics  from  50  Monte  Carlo  runs  of  the  2  state  filters  on 
the  S-tum  target  are  given  in  Figures  C-37  through  C-42.  The  Chi-square  test  of  the  bearing 
innovations  is  given  in  Figure  C-37.  Most  of  the  points  are  within  the  95-percent  confidence 
region,  but  the  time  average  of  these  points  is  approximately  0.9  which  indicates  that  the  qk  is 

slightly  too  large.  In  Figure  C-38,  the  elevation  innovations  are  consistent  until  the  end  when  the 
true  elevation  is  changing  quickly.  The  whiteness  test  of  the  bearing  innovations  (Figure  C-39) 
indicates  the  bearing  innovations  have  a  slight  negative  correlation  before  and  after  the  maneuver 
(which  begins  at  the  10-s  point  and  ends  at  the  15-s  point)  and  have  a  slight  positive  correlation 
during  the  maneuver.  This  test  indicates  that  the  qk  cannot  be  made  smaller  since  the  bearing 
innovations  will  have  a  stronger  correlation  during  the  maneuver.  In  light  of  this  test,  the  reader 
can  appreciate  why  the  qk  was  not  decreased.  The  whiteness  test  for  the  elevation  innovations 

in  Figure  C-40  gives  the  same  indications  that  the  Chi-square  test  of  the  elevation  innovations 
did  in  Figure  C-38.  The  elevation  innovations  are  uncorrelated  until  the  end  of  the  trajectory.  The 
Chi-square  tests  of  the  bearing  and  elevation  estimation  error  are  shown  in  Figures  C-41  and  C- 
42,  respectively.  As  was  the  case  for  both  the  radial  and  the  crossing  targets,  they  are  consistent 
using  a  one-sided  confidence  region.  As  in  the  previous  two  target  cases,  all  of  the  consistency 
tests  passed,  so  the  process  noise  variances  used  for  these  simulations  are  good  values. 
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8.1.4  Jinking  Target  * 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
2  state  filters  on  the  jinking  target  are  shown  in  Figures  C-43  through  C-46.  In  Figures  C-43  and 
C-44,  the  bearing  and  bearing  rate  estimates  follow  the  true  bearing  and  bearing  rate  through  the 
jinking  maneuvers  reasonably  well.  The  qk  was  intentionally  chosen  large  to  allow  significant 

deviations  in  the  true  bearing  similarly  to  the  S-tum  target  discu.  sed  previously.  In  Figure  C-45, 
the  sample  elevation  estimate  looks  similar  to  the  sample  elevation  estimates  for  the  radial  and 
S-tum  targets  shown  in  Figures  C-3  and  C-31,  respectively.  The  jinking  target  elevation  estimate 
in  Figure  C-45  is  more  similar  to  the  radial  target  elevation  estimate  in  Figure  C-3  since  the#* 

for  these  filters  are  nearly  the  same  (from  Table  1,  5x1  O'4  for  the  jinking  target  and  1x10 3  for 
the  radial  target).  Figure  C-46  shows  that  the  elevation  rate  estimate  does  not  have  a  large  lag 
behind  the  true  elevation  rate  as  one  would  except  for  a  filter  with  a  large  qk. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  50  Monte  Carlo 
runs  of  the  2  state  filters  on  the  jinking  target  are  shown  in  Figures  C-47  through  C-50.  In  Figure 
C-47,  the  RMS  bearing  enror  remains  nearly  constant  throughout  the  trajectory  at  about  1.3  mrad 
(a  25-percent  reduction  from  the  1.74  mrad  measurement  noise  standard  deviation).  As  for  the 
S-tum  target,  one  cannot  easily  see  where  the  maneuvers  occur.  In  Figure  C-48,  one  can  see  the 
jinking  maneuvers  toward  the  end  of  the  trajectory  in  the  RMS  bearing  rate  errors.  In  Figure  C- 
49,  the  RMS  elevation  error  is  about  0.9  mrad  until  the  end  of  the  trajectory  which  is  a  reduction 
of  48-percent  from  1.74  mrad.  The  RMS  elevation  rate  error  (Figure  C-50)  is  steady  at  about 
2.2  mrad/s  until  the  end  of  the  trajectory  where  it  grows  quickly. 

Plots  of  the  consistency  test  statistics  from  50  Monte  Carlo  runs  of  the  2  state  filters  on 
the  jinking  target  are  given  in  Figures  C-51  through  C-56.  Tne  Chi-square  test  of  the  bearing 
innovations  is  given  in  Figure  C-51.  Most  of  the  points  are  within  the  95-percent  confidence 
region.  The  time  average  of  the  bearing  innovations  over  the  first  17  s  of  the  trajectory  is  less 
than  the  optimum  value  of  1,  but  the  time  average  from  17  to  22  s  is  very  close  to  1,  which 
indicates  that  the  qk  is  too  large  over  the  first  17  s  and  is  close  to  optimum  from  17  to  22  s.  In 
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Figure  C-52,  the  elevation  innovations  are  consistent  until  the  end  when  the  true  elevation  is 
changing  quickly.  Since  more  points  fall  below  the  confidence  region  than  above,  the  qk  could 

have  been  decreased,  but  this  would  cause  the  elevation  innovations  to  leave  the  confidence 
region  more  quickly  at  the  end  of  the  trajectory.  The  whiteness  test  of  the  bearing  innovations 
(Figure  C-53)  indicates  the  bearing  innovations  have  a  significant  negative  correlation  until  about 
17  s  into  the  trajectory  where  the  magnitude  of  the  bearing  rate  changes  are  becoming  large.  The 
whiteness  test  for  the  elevation  innovations  (Figure  C-54)  shows  they  are  slightly  correlated  until 
the  end  of  the  trajectory  where  they  become  strongly  correlated.  The  Chi-square  tests  of  the 
bearing  and  elevation  estimation  error  are  shown  in  Figures  C-55  and  C-56,  respectively.  As  was 
the  case  for  all  of  the  previous  targets,  they  are  consistent  using  a  one-sided  confidence  region. 
For  the  jinking  target,  the  consistency  tests  passed  (although  not  with  flying  colors),  so  the 
process  noise  variances  used  for  these  simulations  are  good  values. 

8.2  RESULTS  FOR  3  STATE  FILTER 

The  format  of  this  section  is  the  same  as  the  one  used  in  Section  8.1.  Sample  bearing, 
bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the  3  state  filters  are 
presented.  Then,  the  RMSE  from  30  Monte  Carlo  runs  are  presented.  Each  subsection  concludes 
with  plots  of  the  consistency  test  statistics  from  these  simulations  to  verify  that  the  values  chosen 
for  the  process  noise  variances  are  appropriate.  All  of  the  figures  referred  to  in  this  section  are 
located  in  Appendix  D. 

The  results  for  the  bearing  rate  rate  and  elevation  rate  rate  estimates  are  not  presented. 
This  omission  is  justified  for  two  reasons.  First,  this  work  is  intended  to  support  radar-to-ESM 
correlation  studies,  and  most  radar  systems  can  provide  angle  and  angle  rate  estimates  but  not 
angle  rate  rate  estimates.  Secondly,  the  reader  can  interpolate  the  results  presented  for  the  angle 
and  angle  rate  estimates  to  determine  the  approximate  performance  of  the  3  state  filters  for  the 
angle  rate  rate  estimates. 
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8.2.1  Radial  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
3  s.ate  filters  on  the  radial  target  are  shown  in  Figures  D-l,  D-2,  D-3,  and  D-4,  respectively.  The 
bearing  and  bearing  rate  estimates  in  Figures  D-l  and  D-2  are  rather  noisy,  which  illustrates  one 
of  the  disadvantages  to  using  an  excessively  high-order  filter.  The  filter  is  estimating  both  the 
first  and  second  derivatives  of  the  bearing,  and  from  these  figures,  it  is  obvious  that  the  true  first 
and  second  derivatives  are  zero.  Therefore,  the  filter  is  wasting  information  (i.e.,  the 
measurements)  by  trying  to  estimate  states  that  are  zero.  Figures  D-3  and  D-4  show  that  the 
elevation  and  elevation  rate  estimates  are  somewhat  noisy  in  the  first  20  s  of  the  trajectory,  but 
the  looseness  in  the  filter  allows  the  estimates  to  follow  the  quick  increase  of  the  true  elevation 
and  rate  fairly  well. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  30  Monte  Carlo 
runs  of  the  3  state  filters  on  the  radial  target  are  shown  in  Figures  D-5  through  D-8.,  The 
horizontal  lines  in  Figures  D-5  and  D-7  correspond  to  the  standard  deviation  of  the  measurement 
noise  (0.1  deg  or  1.74  mrad)  and  is  provided  as  a  reference  to  determine  the  amount  of  noise 
removed  by  the  filter.  In  Figure  D-5,  the  RMS  bearing  error  has  been  reduced  from  1.74  mrad 
to  about  1.3  mrad,  which  is  a  25-percent  decrease.  In  Figure  D-6,  the  RMS  bearing  rate  error  is 
approximately  8  mrad/s.  In  Figure  D-7,  the  RMS  elevation  error  has  been  reduced  from  1.74 
mrad  to  1.4  mrad,  which  is  a  20-percent  decrease.  In  Figure  D-8,  the  RMS  elevation  rate  error 
remains  at  about  8  mrad/s  until  the  end  of  the  trajectory  where  the  sharp  rise  in  the  true  elevation 
rate  causes  the  errors  to  increase  quickly. 

Plots  of  the  consistency  test  statistics  from  30  Monte  Carlo  runs  of  the  3  state  filters  on 
the  radial  target  are  given  in  Figures  D-9  through  D-14.  The  area  between  the  horizontal  lines 
is  the  95-percent  confidence  region  (refer  to  Section  8  1  for  an  explanation  of  the  Y  axis  labels 
used  in  these  figures).  The  Chi-square  test  of  the  bearing  innovations  is  given  in  Figure  D-9. 
Most  of  the  points  are  within  the  95-percent  confidence  region,  so  the  bearing  innovations  are 
consistent  with  their  covariances.  Since  the  time  average  of  these  points  is  close  to  0.8,  th aqk 
could  actually  be  slightly  smaller  than  it  is.  The  Chi-square  test  of  the  elevation  innovations 
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presented  in  Figure  D-10  shows  that  the  elevation  innovations  are  consistent  The  whiteness  test 
of  the  bearing  innovations  is  shown  in  Figure  D-ll.  For  this  figure  and  all  subsequent  whiteness 
tests,  only  the  for  k-j=l  is  shown.  As  shown  in  this  figure,  the  innovations  have  a  slight 

negative  correlation,  which  is  another  indication  that  the  qk  could  have  been  smaller.  The 

whiteness  test  for  the  elevation  innovations  shown  in  Figure  D-12  is  similar  to  that  of  the  bearing 
innovations.  They  have  a  slight  negative  correlation.  The  Chi-square  test  of  the  bearing  estimation 
error  is  shown  in  Figure  D-13.  The  confidence  region  is  not  shown  since  the  thresholds  are  on 
the  order  of  one.  As  was  the  case  for  the  2  state  filters,  the  3  state  filters  are  usually  pessimistic. 
For  this  reason,  a  one-sided  confidence  region  will  be  used  on  most  of  the  following  Chi-square 
tests  of  the  state  estimation  error.  The  Chi-square  test  of  the  elevation  estimation  error  is  shown 
in  Figure  D-14.  The  data  up  to  24  s  through  the  trajectory  are  on  the  order  of  10'3  and  cannot  be 
seen  in  the  figure.  The  large  errors  in  the  elevation  rate  rate  estimates  are  the  major  cause  for  the 
sharp  increase  in  ek  in  the  last  second  of  Figure  D-14.  Again,  the  missile  is  at  a  very  short  range 

at  this  point,  so  large  elevation  estimate  errors  are  not  significant  Thus,  all  of  the  consistency 
tests  passed,  so  the  process  noise  variances  used  for  these  simulations  are  good  values. 

8.2.2  Crossing  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
3  state  filters  on  the  crossing  target  are  shown  in  Figures  D-15  through  D-18,  respectively.  In 
Figure  D-15,  the  difference  between  the  sample  bearing  estimate  and  the  true  bearing  is  too  small 
to  be  seen.  The  errors  in  the  sample  bearing  rate  estimate  in  Figure  D-16  are  about  the  same 
magnitude  through  the  entire  trajectory.  The  sample  elevation  and  elevation  rate  estimates  shown 
in  Figures  D-17  and  D-18  are  somewhat  noisy  and  are  similar  to  the  bearing  and  bearing  rate 
estimates  for  the  radial  target  shown  in  Figures  D-3  and  D-4. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  30  Monte  Carlo 
runs  of  the  3  state  filters  on  the  crossing  target  are  shown  in  Figures  D-19,  D-20,  D-21,  and  D- 
22,  respectively.  In  Figure  D-19,  the  RMS  bearing  error  is  approximately  0.9  mrad  throughout 
the  trajectory,  which  is  a  48-percent  reduction  from  the  1.74  mrad  measurement  noise  standard 


27 


NSWCDD/TR-92/503 


deviation.  The  RMS  bearing  rate  error  remains  virtually  constant  at  1.8  mrad/s.  There  is  a  slight 
increase  in  the  RMS  error  around  the  20-s  point,  which  is  a  result  of  the  true  bearing  rate  rate 
approaching  zero  at  this  point.  In  Figure  D-21,  the  RMS  elevation  error  has  been  reduced  from 
1.74  mrad  to  0.8  mrad  which  is  a  54-percent  reduction.  The  RMS  elevation  rate  error  remains 
at  1.5  mrad/s  throughout  the  trajectory  as  shown  in  Figure  D-22. 

Plots  of  the  consistency  test  statistics  from  30  Monte  Carlo  runs  of  the  3  state  filters  on 
the  crossing  target  are  given  in  Figures  D-23  through  D-28.  All  of  these  consistency  tests  passed, 
so  the  chosen  values  for  qk  are  appropriate. 

8.2.3  S-tum  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
3  state  filters  on  the  S-tum  target  are  shown  in  Figures  D-29  through  D-32.  In  Figure  D-30,  the 
largest  error  in  the  sample  bearing  rate  estimate  occurs  at  about  the  13-s  point  where  the  bearing 
rate  rate  changes  quickly. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  30  Monte  Carlo 
runs  of  the  3  state  filters  on  the  S-tum  target  are  shown  in  Figures  D-33  thrcugn  D-36.  In  Figure 
D-33,  the  RMS  bearing  error  remains  nearly  constant  throughout  the  trajectory  at  about  1.1  mrad 
(a  37-percent  reduction  from  the  1.74  mrad  measurement  noise  standard  deviation).  The 
maneuver  that  occurs  near  the  middle  of  the  trajectory  cannot  be  seen  in  the  RMS  bearing  errors. 
In  Figure  D-34,  the  maneuver  is  easier  to  see  in  the  RMS  bearing  rate  errors  which  vary 
between  3  mrad/s  and  6  mrad/s.  In  Figure  D-35,  the  RMS  elevation  error  is  about  0.8  mrad 
throughout  the  trajectory,  which  is  a  reduction  of  54-percent  from  1.74  mrad.  At  the  end  of  the 
trajectory,  the  elevation  errors  are  growing  but  do  not  exceed  the  1.74-mrad  line.  The  RMS 
elevation  rate  error  (Figure  D-36)  is  steady  at  1  mrad/s  until  the  end  where  it  begins  increasing. 

Plots  of  the  consistency  test  statistics  from  30  Monte  Carlo  runs  of  the  3  state  filters  on 
the  S-tum  target  are  given  in  Figures  D-37  through  D-42.  The  Chi-square  tests  of  the  bearing  and 
elevation  estimation  erroi  are  shown  in  Figures  D-41  and  D-42,  respectively.  As  was  the  case  for 
both  the  radial  and  the  crossing  targets,  they  are  consistent  using  a  one-sided  confidence  region. 
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The  large  peaks  liiat  appear  in  these  figures  are  primarily  due  to  the  large  errors  in  the  angle  rate 
rates.  As  in  the  previous  two  target  cases,  all  of  the  consistency  tests  passed,  so  the  process  noise 
variances  used  for  these  simulations  are  good  values. 

8.2.4  Jinking  Target 

The  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  a  single  run  of  the 
3  state  filters  on  the  jinking  target  are  shown  in  Figures  D-43  through  D-46.  In  Figures  D-43  and 
D-44,  the  bearing  and  bearing  rate  estimates  follow  the  true  bearing  and  bearing  rate  through  the 
jinking  maneuvers  reasonably  well.  In  Figure  D-45,  the  sample  elevation  estimate  looks  similar 
to  the  sample  elevation  estimates  for  the  raaial  and  S-tum  targets  shown  in  Figures  D-3  and  D- 
31,  respectively. 

The  RMS  bearing,  bearing  rate,  elevation,  and  elevation  rate  errors  from  30  Monte  Carlo 
runs  of  the  3  state  filters  on  the  jinking  target  are  shown  in  Figures  D-47  through  D-50.  In  Figure 
D-47,  the  RMS  bearing  error  is  approximately  1.3  mrad  (a  25-percent  reduction  from  the  1.74 
mrad  measurement  noise  standard  deviation).  In  Figure  D-48,  one  can  see  the  jinking  maneuvers 
toward  the  end  of  the  trajectory  in  the  RMS  bearing  rate  errors.  In  Figure  D-49,  the  RMS 
elevation  error  is  about  1  mrad  until  the  end  of  the  trajectory  which  is  i  eduction  of  43-percent 
from  1.74  mrad.  The  RMS  elevation  rate  error  (Figure  D-50)  is  steady  at  about  2.2  mrad/s  until 
the  end  of  the  trajectory  where  it  grows  quickly.. 

Plots  of  the  consistency  test  statistics  from  30  Monte  Carlo  runs  of  the  3  state  fibers  on 
the  jinking  target  are  given  in  Figures  D-51  through  D-56.  All  the  consistency  tests  passed,  so 
the  process  noise  variances  used  for  these  simulations  are  good  values. 

8.3  COMPARISON  OF  2  AND  3  STATE  FILTERS 

The  reader  may  have  noticed  that  little  effort  has  been  made  to  compare  the  results  of  the 
2  and  3  state  filters  in  the  previous  sections.  This  section  was  included  to  fill  this  void.  Monte 
Carlo  simulations  were  performed  and  the  statistic  given  by  Equation  (4.1)  was  computed  to 
compare  the  bearing,  bearing  rate,  elevation,  and  elevation  rate  estimates  from  the  2  and  3  state 
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filters.  The  results  of  these  simulations  are  presented  in  this  section.  All  of  the  figures  referred 
to  in  this  section  are  located  in  Appendix  E. 

In  all  of  the  figures,  if  the  comparison  statistic  is  above  the  upper  threshold,  then  the 
hypothesis  that  the  2  state  filter’s  estimates  are  better  than  the  3  state  filter’s  estimates  is 
accepted.  Likewise,  if  the  comparison  statistic  is  below  the  lower  tnreshold,  then  the  hypothesis 
that  the  3  state  filter’s  estimates  are  better  than  the  2  state  filter’s  estimates  is  accepted.  If  the 
comparison  statistic  is  between  the  two  thresholds,  then  neither  hypothesis  is  accepted.  Both 
thresholds  correspond  to  95-percent  probabilities  that  the  hypothesis  is  correct. 

8.3.1  Radial  Target 

Plots  of  the  results  of  a  comparison  between  the  2  and  3  state  filters  for  50  Monte  Carlo 
runs  using  the  radial  target  are  shown  in  Figures  E-l  to  E-4.  The  performance  of  the  2  state  filter 
was  significantly  better  than  that  of  the  3  state  filter  since  all  of  the  test  statistics  are  greater  than 
the  upper  threshold  that  defines  the  95-percent  confidence  region.  The  only  exception  is  the 
elevation  statistics  in  the  last  few  seconds  of  the  trajectory  (Figures  E-3  and  E-4).  Here  the  true 
elevation  has  a  significant  second  derivative,  and  the  3  state  filter’s  estimates  are  better. 

8.3.2  Crossing  Target 

Plots  of  the  results  of  a  comparison  between  the  2  and  3  state  filters  for  50  Monte  Carlo 
runs  using  the  crossing  target  are  shown  in  Figures  E-5  to  E-8.  In  Figures  E-5  and  E-6,  the 
bearing  and  bearing  rate  comparison  statistics  are  below  the  lower  threshold  for  most  of  the 
trajectory,  so  the  3  state  filter’s  bearing  and  bearing  rate  estimates  are  better.  Notice  at  the  21-s 
point  of  the  trajectory  that  the  2  state  filter’s  estimates  are  better  since  the  true  bearing  rate  rate 
is  approximately  zero.  From  Figures  E-7  and  E-8,  the  2  state  filter’s  elevation  and  elevation  rate 
estimates  are  better  throughout  the  trajectory. 

8.3.3  S-tum  Target 

Plots  of  the  results  of  a  comparison  between  the  2  and  3  state  filters  for  50  Monte  Carlo 
runs  using  the  S-tum  target  are  shown  in  Figures  E-9  to  E-l 2.  Figures  E-9  and  E-xO  show  that 
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both  filters  give  bearing  and  bearing  rate  estimates  of  approximately  the  same  quality.  In  the  first 
turn  on  the  S-tum,  which  occurs  from  10  s  to  12.5  s  in  these  figures,  the  3  state  filter’s  estimates 
appear  to  e  better.  However,  not  enough  points  fall  below  the  lower  threshold  to  justify  a 
conclusion  that  the  3  state  filter’s  estimates  are  superior.  In  Figures  E-ll  and  E-12,  all  of  the 
points  before  the  20  s  point  are  above  the  upper  threshold,  so  the  2  state  filter’s  estimates  are 
better.  The  3  state  filter’s  estimates  are  definitely  better  during  the  last  few  seconds  of  the 
trajectory. 

8.3.4  Jinking  Target 

Plots  of  the  results  of  a  amparison  between  the  2  and  3  state  filters  for  50  Monte  Carlo 
runs  using  the  jinking  target  are  shown  in  Figures  E-13  to  E-16.  In  Figures  E-13  and  E-14,  most 
of  the  poinu'  are  above  the  upper  threshold  until  the  last  few  maneuvers  occur  from  17s  to  21 
s.  Thus,  the  2  state  filter’s  bearing  and  bearing  rate  estimates  are  bettcr  until  the  end  of  the 
trajectory.  In  Figures  E-15  and  E-16,  most  of  the  points  fall  between  the  two  thresholds,  so 
neither  filter’s  elevation  and  elevation  rate  estimates  are  better.  Of  course,  in  the  last  few  seconds 
of  the  trajectory,  the  3  state  filter’s  estimates  are  better. 
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CHAPTER  9 

CONCLUSIONS  AND  RECOMMENDATIONS 

Two  Kalman  filters  were  tested  for  estimating  the  true  bearing,  nearing  rate,  elevation, 
and  elevation  rate  of  high-speed  targets  using  simulated  measurements  from  a  shipboard  ESM 
receiver.  The  two  Kalman  filters  tested  were  described  in  terms  of  the  state  transition  matrix.  The 

2  state  filter  predicts  the  true  state  using  a  constant  first  derivative  assumption,  and  the  3  state 
filter  predicts  the  true  state  using  a  constant  second  derivative  assumption.  The  steady  state 
versions  of  these  filters  are  the  well  known  a,Pand  an  a, p,y  filters.  The  steady  state  filters 
would  be  used  in  an  actual  implementation  since  they  are  computationally  simpler.  The  Kalman 
filter  versions  were  used  in  the  research  to  simplify  the  coding  of  the  algorithms. 

Consistency  tests  were  described  as  an  effective  method  to  determine  the  best  process 
noise  variance  (i.e.,  the  steady  state  gains)  for  each  filter  for  each  target  trajectory  tested.  An 
example  of  this  method  of  process  noise  variance  selection  was  presented. 

The  four  target  trajectories  used  to  test  the  tracking  filters  were  presented.  All  were  Mac.1. 

3  targets  with  a  data  rate  of  10  Hz  and  a  measurement  noise  standard  deviation  of  0.1  deg  (1.74 
mrad).  Two  targets  moved  in  straight  lines,  and  the  other  two  performed  maneuvers  with  15-g 
accelerations. 

A  statistic  that  can  be  used  to  compare  the  estimates  generated  by  two  tracking  filters  was 
presented.  This  statistic  was  used  to  compare  the  bearing,  bearing  rate,  elevation,  and  elevation 
rate  estimates  from  the  2  and  3  state  filters  for  all  four  target  trajectories.  The  2  state  filter  gave 
better  bearing  and  bearing  rate  estimates  for  all  targets  except  the  crossing  target.  The  3  state 
filter’s  bearing  and  bearing  rate  estimates  were  superior  for  most  of  the  crossing  trajectory.  The 
2  state  filter  also  gave  better  elevation  and  elevation  rate  estimates  for  all  of  the  targets. 
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However,  for  the  radial,  S-tum,  and  jinking  targets,  the  3  state  filter’s  elevation  and  elevation  rate 
estimates  were  better  in  the  last  few  seconds  of  the  trajectories  since  the  true  elevation  rate  was 
changing  quickly. 

Both  filters  are  recommended  for  radar-to-ESM  track  correlation  research.  The  2  state 
filter  should  be  tested  first  since  it  is  simpler  and  provides  better  estimates  for  most  of  the  cases 
tested. 
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APPENDIX  A 

PLOTS  OF  TARGET  TRAJECTORIES 

This  appendix  presents  the  four  target  trajectories  used  to  test  the  2  and  3  state  filters.  The 
trajectories  are  presented  in  the  following  order:  radial  target,  crossing  target,  S-tum  target,  and 
jinking  target. 
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FIGURE  A-l.  TRUE  BEARING  FOR  RADIAL  TARGET 
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FIGURE  A-2.  TRUE  ELEVATION  FOR  RADIAL  TARGET 
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FIGURE  A  ll,  TRUE  ACCELERATION  FOR  S-TURN  TARGET 
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FIGURE  A- 12.  TRUE  BEARING  FOR  S-TURN  TARGET 
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FIGURE  A- 13.  TRUE  BEARING  RATE  FOR  S-TURN  TARGET 


FIGURE  A- 14.  TRUE  ELEVATION  FOR  S-TURN  TARGET 
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FIGURE  A-15.  TRUE  ELEVATION  RATE  FOR  S-TURN  TARGET 
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FIGURE  A- 19.  TRUE  BEARING  FOR  JINKING  TARGET 
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FIGURE  A-22.  TRUE  ELEVATION  RATE  FOR  JINKING  TARGET 
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APPENDIX  B 

PLOTS  OF  PROCESS  NOISE  VARIANCE  SELECTION 

This  appendix  presents  the  results  of  the  example  of  the  process  noise  variance  selection 
procedure.  This  example  is  for  the  3  state  filter  on  the  crossing  target.  The  consistency  tests  are 
used  to  determine  the  best  value  of  the  process  noise  variance  for  the  particular  trajectory. 
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FIGURE  B-I.  CHI-SQUARE  TEST  OF  BEARING  INNOVATIONS  FOR  3 
STATE  FILTER  ON  CROSSING  TARGET  WITH  q=106 


FIGURE  B-2.  WHITENESS  TEST  OF  BEARING  INNOVATIONS  FOR  3  STATE 
FILTER  ON  CROSSING  TARGET  WITH  q=ltf6 
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FIGURE  B-3.  CHI-SQUARE  TEST  OF  BEARING  INNOVATIONS  FOR  3 
STATE  FILTER  ON  CROSSING  TARGET  WITH  q=107 


FIGURE  B-4,  WHITENESS  TEST  OF  BEARING  INNOVATIONS  FOR  3  STATE 
FIL1ER  ON  CROSSING  TARGET  WITH  q=10'7 
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APPENDIX  C 

PLOTS  OF  2  STATE  FILTER  RESULTS 

This  appendix  presents  the  results  of  simulations  performed  for  the  2  state  filters.  The 
results  are  presented  in  the  following  order:  radial  target,  crossing  target,  S-turn  target,  and 
jinking  target. 
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FIGURE  C-2.  SAMPLE  BEARING  RATE  ESTIMATE  FOR  2  STATE  FILTER 

ON  RADIAL  TARGET 
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1GURE  C-3, 


SAMPLE  ELEVATION  ESTIMATE  FOR  2  STATE  FILTER  ON 
RADIAL  TARGET 
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FIGURE  C-4.  SAMPLE  ELEVATION  RATE  ESTIMATE  FOR  2  STATE  FILTER 

ON  RADIAL  TARGET 
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FIGURE  C-8.  RMS  ELEVATION  RATE  ERROR  FOR  2  STATE  FILTER  ON 

RADIAL  TARGET 
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FIGURE  C-10.  CHI-SQUARE  TEST  OF  ELEVATION  INNOVATIONS  FOR  2 
STATE  FILTER  ON  RADIAL  TARGET 
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FIGURE  C-17.  SAMPLE  ELEVATION  ESTIMATE  FOR  2  STATE  FILTER  ON 

CROSSING  TARGET 
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FIGURE  C-18,  SAMPLE  ELEVATION  RATE  ESTIMATE  FOR  2  STATE 
FILTER  ON  CROSSING  TARGET 
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FIGURE  C-20.  RMS  BEARING  RATE  ERROR  FOR  2  STATE  HLTER  ON 

CROSSING  TARGET 
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FIGURE  C-22.  RMS  ELEVATION  RATE  ERROR  FOR  2  STATE  FILTER  ON 

CROSSING  TARGET 
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FIGURE  C-25.  WHITENESS  TEST  OF  BEARING  INNOVATIONS  FOR  2 
STATE  FILTER  ON  CROSSING  TARGET 


FIGURE  C-26.  WHITENESS  TEST  OF  ELEVATION  INNOVATIONS  FOR  2 
STATE  FILTER  ON  CROSSING  TARGET 
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FIGURE  C-31.;  SAMPLE  ELEVATION  ESTIMATE  FOR  2  STATE  FILTER  ON 

S-TURN  TARGET 
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FIGURE  C-32.  SAMPLE  ELEVATION  RATE  ESTIMATE  FOR  2  STATE 
FILTER  ON  S-TURN  TARGET 
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FIGURE  C-35.  RMS  ELEVATION  ERROR  FOR  2  STATE  FILTER  ON  S-TURN 

TARGET 


FIGURE  C-36.  RMS  ELEVATION  RATE  ERROR  FOR  2  :>  i  VIE  FILTER  ON  S- 

TURN  TARGET 
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FIGURE  C-38.  CHI-SQUARE  TEST  OF  ELEVATION  INNOVATIONS  FOR  2 
STATE  FILTER  ON  S-TURN  TARGET 
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[GURE  C41,  CHI-SQUARE  TEST  OF  BEARING  ESTIMATION  ERROR  FOR 
2  STATE  FILTER  ON  S-TURN  TARGET 
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FIGURE  CHI-SQUARE  TEST  OF  ELEVATION  ESTIMATION  ERROR 
FOR  2  STATE  FILTER  ON  S-TURN  TARGET 
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FIGURE  C-43.  SAMPLE  BEARING  ESTIMATE  FOR  2  STATE  FILTER  ON 

JINKING  TARGET 
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FIGURE  C-45,  SAMPLE  ELEVATION  ESTIMATE  FOR  2  STATE  FILTER  ON 

JINKING  TARGET 
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FIGURE  C46.  SAMPLE  ELEVATION  RATE  ESTIMATE  FOR  2  STATE 
FILTER  ON  JINKING  TARGET 
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FIGURE  C-50.  RMS  ELEVATION  RATE  ERROR  FOR  2  STATE  FILTER  ON 

JINKING  TARGET 
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FIGURE  C-56,  CHI-SQUARE  TEST  OF  ELEVATION  ESTIMATION  ERROR 
FOR  2  STATE  FILTER  ON  JINKING  T  ARGET 
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APPENDIX  D 

PLOTS  OF  3  STATE  FILTER  RESULTS 

This  appendix  presents  the  results  of  simulations  performed  fo,  .he  3  state  filters.  The 
results  are  presented  in  the  following  order:  radial  target,  crossing  target,  S-tum  target,  and 
jinking  target. 
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FIGURE  D-3.  SAMPLE  ELEVATION  ESTIMATE  FOR  3  STATE  FILTER  ON 

RADIAL  TARGET 
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ON  RADIAL  TARGET 
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FIGURE  D-ll.  WHITENESS  TEST  OF  BEARING  INNOVATIONS  FOR  3 
STATE  FILTER  ON  RADIAL  TARGET 


FIGURE  D-12.  WHITENESS  TEST  OF  ELEVATION  INNOVATIONS  FOR  3 
STATE  FILTER  ON  RADIAL  TARGET 
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FIGURE  D-14.  CHI-SQUARE  TEST  OF  ELEVATION  ESTIMATION  ERROR 
FOR  3  STATE  FILTER  ON  RADIAL  TARGET 
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FIGURE  D-15.  SAMPLE  BEARING  ESTIMATE  FOR  3  STATE  FILTER  ON 

CROSSING  TARGET 


FIGURE  D-16.  SAMPLE  BEARING  RATE  ESTIMATE  FOR  3  STATE  FILTER 

ON  CROSSING  TARGET 
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IGURE  D-17..  SAMPLE  ELEVATION  ESTIMATE  FOR  3  STATE  FILTER  0. 

CROSSING  TARGET 
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FIGURE  D-18.  SAMPLE  ELEVATION  RATE  ESTIMATE  FOR  3  STATE 
FILTER  ON  CROSSING  TARGET 
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FIGURE  D-20.  RMS  BE/JRING  RATE  ERROR  FOR  3  STATE  FILTER  ON 

CROSSING  TARGET 
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FIGURE  D-23. 


CHI-SQUARE  TEST  OF  BEARINC 
STATE  FILTER  ON  CROSSING  T 
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FIGURE  D-24. 


CHI-SQUARE  TEST  OF  ELEVATIC 
STATE  FILTER  ON  CROSSING  T 
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FIGURE  D-25.  WHITENESS  TEST  OF  BEARING  INNOVATIONS  FOR  3 
STATE  FILTER  ON  CROSSING  TARGET 


FIGURE  D-26.  WHITENESS  TEST  OF  ELEVATION  INNOVATIONS  FOR 
STATE  FILTER  ON  CROSSING  TARGET 
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FIGURE  D-30.  SAMPLE  BEARING  RATE  ESTIMATE  FOR  3  STATE  FILTER 

ON  S-TURN  TARGET 
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FIGURE  D-34.  RMS  BEARING  RATE  ERROR  FOR  3  STATE  FILTER  ON  S- 

TURN  TARGET 
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FIGURE  D-38.  CHI-SQUARE  TEST  OF  ELEVATION  INNOVATIONS  FOR  3 
STATE  FILTER  ON  S-TURN  TARGET 
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FIGURE  D-39.  WHITENESS  TEST  OF  BEARING  INNOVATIONS  FOR  3 
STATE  FILTER  ON  S-TURN  TARGET 


FIGURE  D-40.  WHITENESS  TEST  OF  ELEVATION  INNOVATIONS  FOR  3 
STATE  FILTER  ON  S-TURN  TARGET 
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FIGURE  D-41.  CHI-SQUARE  TEST  OF  BEARING  ESTIMATION  ERROR  FOR 
3  STATE  FILTER  ON  S-TURN  TARGET 
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1GURE  D-45.  SAMPLE  ELEVATION  ESTIMATE  FOR  3  STATE  FILTER  ON 

JINKING  TARGET 
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FIGURE  D46.  SAMPLE  ELEVATION  RATE  ESTIMATE  FOR  3  STATE 
FILTER  ON  JINKING  TARGET 
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FIGURE  D-48.  RMS  BEARING  RATE  ERROR  FOR  3  STATE  FILTER  ON 

JINKING  TARGET 
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FIGURE  D-50.  RMS  ELEVATION  RATE  ERROR  FOR  3  STATE  FILTER  ON 

JINKING  TARGET 
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FIGURE  D-55.  CHI-SQUARE  TEST  OF  BEARING  ESTIMATION  ERROR  FOR 
3  STATE  FILTER  ON  JINKING  TARGET 
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FIGURE  D-56.  CHI-SQUARE  TEST  OF  ELEVATION  ESTIMATION  ERROR 
FOR  3  STATE  FILTER  ON  JINKING  TARGET 
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APPENDIX  E 

PLOTS  OF  COMPARISON  OF  2  AND  3  STATE  FILTERS 

This  appendix  presents  the  results  of  simulations  performed  to  compare  the  performance 
of  the  2  and  3  state  filters.  The  results  are  presented  in  the  following  order:  radial  target,  crossing 
target,  S-tum  target,  and  jinking  target 
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FIGURE  E-l..  COMPARISON  OF  BEARING  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  RADIAL  TARGET 


FIGURE  E-2.  COMPARISON  OF  BEARING  RATE  ESTIMATES  FROM  2  AND 
3  STATE  FILTERS  ON  RADIAL  TARGET 
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FIGURE  E-3.  COMPARISON  OF  ELEVATION  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  RADIAL  TARGET 
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FIGURE  E-4.  COMPARISON  OF  ELEVATION  RATE  ESTIMATES  FROM  2 
AND  3  STATE  FILTERS  ON  RADIAL  TARGET 
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FIGURE  E-5,  COMPARISON  OF  BEARING  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  CROSSING  TARGET 
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FIGURE  E-6.  COMPARISON  OF  BEARING  RATE  ESTIMATES  FROM  2  AND 
3  STATE  FILTERS  ON  CROSSING  TARGET 
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FIGURE  E-7.  COMPARISON  OF  ELEVATION  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  CROSSING  TARGET 


FIGURE  E-8.  COMPARISON  OF  ELEVATION  RATE  ESTIMATES  FROM  2 
AND  3  STATE  FILTERS  ON  CROSSING  TARGET 
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FIGURE  E-9.  COMPARISON  OF  BEARING  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  S-TURN  TARGET 


0  5  10  15  20  25 


TIME  (S) 

FIGURE  E-10.  COMPARISON  OF  BEARING  RATE  ESTIMATES  FROM  2 
AND  3  STATE  FILTERS  ON  S-TURN  TARGET 
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FIGURE  E-ll.  COMPARISON  OF  ELEVATION  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  S-TURN  TARGET 


FIGURE  E-12.  COMPARISON  OF  ELEVATION  RATE  ESTIMATES  FROM  2 
AND  3  STATE  FILTERS  ON  S-TURN  TARGET 
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FIGURE  E-13,  COMPARISON  OF  BEARING  ESTIMATES  FROM  2  AND  3 
STATE  FILTERS  ON  JINKING  TARGET 
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FIGURE  E-14.  COMPARISON  OF  BEARING  RATE  ESTIMATES  FROM  2 
AND  3  STATE  FILTERS  ON  JINKING  TARGET 
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1GURE  E-15.  COMPARISON  OF  ELEVATION  ESTIMATES  FROM  2  AND 
STATE  FILTERS  ON  JINKING  TARGET 


FIGURE  E-16.  COMPARISON  OF  ELEVATION  RATE  ESTIMATES  FROM  2 
AND  3  STATE  FILTERS  ON  JINKING  TARGET 
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